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Introduction 


The  objective  of  the  study  is  to  develop  a  multimodal  approach  to  assess  concussion  in  the 
acute  phase  of  injury.  Several  different  measures  are  included  in  the  assessment.  These  include 
functional  brain  measurement  techniques,  anatomical  brain  mapping  techniques,  biomarkers 
(blood  protein  levels),  neurocognitive  performance  (computer-based  tests),  balance  testing, 
and  behavioral  health  measures.  The  ultimate  aim  is  to  develop  an  objective  method  to 
diagnose  mTBI  in  military  personnel  proximal  to  the  time  of  injury  and  facilitate  'return  to  duty' 
decisions. 

Keywords 

Traumatic  brain  injury,  mTBI,  concussion,  magnetoencephalography,  MEG,  MRI,  blood 
biomarkers,  actigraphy 

Accomplishments 

Major  Goals 

Determine  the  pattern  of  physiological  and  behavioral  markers  that  characterize  mild  traumatic 
brain  injury  (mTBI)  to  develop  a  comprehensive  methodology  for  the  diagnosis  and  assessment 
of  this  disorder. 

Milestones  were  to  consent  50  mTBI  patients  and  50  non-head  trauma  controls  per  year,  for  3 
multimodal  studies  within  the  first  month  of  injury. 

Accomplishments 

Recruitment  in  our  research  project  finally  started  in  September  2014.  Interactions  with  the 
emergency  department  medical  staff  have  developed  successfully  based  on  frequent  visits  from 
HMRI  to  raise  awareness.  Weekly  referrals  have  been  delayed  during  times  of  our  laboratory 
down  time,  due  to  instrument  repairs  and  personnel  absences).  We  were  referred  a  total  of  25 
patients  from  the  Hospital,  and  consented  a  total  of  15  patients,  3  of  whom  have  dropped  out. 
Six  with  mTBI  and  six  non-head  trauma  controls  have  been  studied.  This  patient  number  that 
have  been  compliant  with  the  full  study  represents  accomplishment  of  25%  of  our  milestones 
for  the  period  thus  far. 

The  first  analysis  compared  behavioral  and  physiological  data  between  acute  mTBI  patients  and 
non-head  trauma  controls  in  the  first  days  after  injury.  Candidate  blood  biomarker  methods 
have  been  developed  with  femtomole  sensitivity  for  multiple  peptides  for  nine  proteins  that 
may  arise  from  brain  trauma,  using  the  newly  acquired  mass  spectrometer;  analyses  are  in 
progress.  A  larger  sample  size  is  required  before  their  analysis.  Multiple  modalities  of 
behavioral,  electrophysiological,  and  most  strikingly,  MEG  changes  have  been  found  in  this 
preliminary  study.  The  MEG  of  all  (n  =  6)  mTBI  patients  had  delta  activity  in  the  frontal  lobes 
that  was  absent  in  all  controls  (n  =  6).  A  scientific  abstract  on  these  findings  has  been  accepted 
for  publication  and  these  results  may  be  the  forerunner  of  an  objective  test  for  acute 
concussion.  More  complete  details  on  the  methods,  analysis,  results,  and  images  are  in  the 
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appendix  abstract  included  with  this  report.  Analysis  of  this  same  patient  group  one  month 
after  injury  when  symptoms  after  mTBI  had  resolved,  is  ongoing.  This  one-month  data  may 
provide  evidence  of  neuronal  damage  recovery  from  the  trauma;  objective  demonstration  of 
improvement  would  have  exciting  potential  for  predicting  and  monitoring  recovery  after  acute 
brain  trauma,  compared  to  the  persistent  delta  waves  reported  in  the  majority  of  chronic  TBI 
patients. 

Training  and  Professional  Development 

Because  of  the  personnel  changes,  four  key  personnel  acquired  critical  training  as  follows: 

•  Dr.  Xianghong  Arakaki  MD,  PhD,  is  an  experienced  electrophysiologist  but  she  had  not 
studied  human  research  previously.  This  year,  she  has  been  trained  by  multiple  local  and 
national  experts  and  has  become  exceedingly  proficient  at  resting  EEG,  auditory  and  executive 
function  evoked  potentials,  and  MEG  analysis.  This  has  profoundly  changed  her  career 
direction,  in  a  positive  way. 

•  Ms.  Thao  Tran  is  an  expert  MR  spectroscopy  and  imaging  technician  and  gained  expertise  at 
UCSD  and  UCSF  to  acquire  quality  MEG  data.  She  has  also  trained  an  HMRI  technician,  David 
Strickland,  to  maintain  the  MEG  instrument  and  perform  routine  MEG  data  acquisition. 

•  Cherise  Charleswell,  MPH,  had  been  trained  by  our  local  neuropsychology  team  to 
administer  all  of  the  neuropsychological  and  behavioral  tests. 

•  The  PI  required  and  received  intensive  training  to  administer  this  project,  learn  the 
additional  components  of  the  science  of  the  project.  This  was  possible  only  with  the  extensive 
support  from  HMRI,  USCD,  UCSF,  and  University  of  Houston.  Dr.  Harrington  was  also  required 
to  reduce  activities  in  other  research  duties  to  devote  additional  time  to  this  mTBI  project. 

Results  Dissemination 

An  abstracted  titled  "Brain  Activation  Profiles  in  mTBI:  Evidence  from  Combined  Resting-State 
EEG  and  MEG  Activity"  has  been  accepted  for  oral  presentation  at  the  37th  Annual 
International  Conference  of  the  IEEE  Engineering  in  Medicine  and  Biology  Society  to  be  held  in 
MiCo,  Milano  Conference  Center,  Milano,  Italy  during  August  25-29,  2015.  Authors  Lianyang  Li, 
Mattia  Federico  Pagnotta,  Xianghong  Arakaki,  Thao  Tran,  David  Strickland,  Michael  Harrington, 
George  Zouridakis 

Plans  to  Accomplish  Goals  in  Next  Reporting  Period 

We  plan  to  actively  continue  to  recruit  mTBI  patients  and  controls  and  extend  the  analysis.  The 
highest  priority  is  to  compare  the  clinical  improvement  of  mTBI  patients  with  objective 
biomarkers;  in  particular,  we  plan  to  closely  analyze  the  MEG  delta  abnormalities  found  in  all 
mTBI  patients  at  the  30-day  period.  We  are  considering  reassessing  these  patients  and  controls 
at  an  additional  time  point  such  as  at  six  months  to  determine  how  these  MEG  observations 
change.  We  will  seek  opinion  from  our  army  advisors  before  making  this  decision  once  we  have 
completed  the  initial  30-day  analysis. 

Impact 

Nothing  to  report. 
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Changes/Problems 


Major  personnel  changes  in  the  summer  of  2014  included: 

•  Captain  Merrifield  left  this  project  in  July  2014  and  his  major  role  has  required  considerable 
adjustments.  The  unexpected  requirements  for  acquiring  at  least  some  of  his  neurobehavioral 
and  MEG  expertise  were  challenging  but  have  been  accomplished. 

•  Dr.  Brian  Ross  retired  in  July  2014.  He  has  been  a  world  leading  in  MR  technologies.  The 
essential  MR  component  for  this  project  has  been  maintained  fully  by  Thao  Tran,  trained  by  Dr. 
Ross. 

•  Analysis  of  EEG,  MEG,  and  MRI  data  with  source  localization  and  co-registration  has  been 
achieved  by  a  combination  of  Dr.  Arakaki,  Ms.  Thao  Tran,  and  our  external  consultant.  Dr. 
George  Zouridakis  in  Houston,  Texas. 

Products 

Nothing  to  report. 

Participants  &  Other  Collaborating  Organizations 


Name 

George  Zouridakis 

Project  Role 

Consultant  Co-Investigator 

Researcher  Identifier  (e.g.  ORCID  ID): 

0000-0002-7770-9857 

Nearest  person  month  worked 

2 

Contribution  to  Project 

Analysis  of  EEG,  MEG,  and  MRI  data  with 
source  localization  and  co-registration 

Funding  Support 

Subcontract  from  HMRI 

Name 

Lianyang  Li 

Project  Role 

Graduate  Student 

Researcher  Identifier  (e.g.  ORCID  ID): 

0000-0002-9600-3120 

Nearest  person  month  worked 

6 

Contribution  to  Project 

Analysis  of  EEG,  MEG,  and  MRI  data  with 
source  localization  and  co-registration 

Funding  Support 

Subcontract  from  HMRI 

Special  Reporting  Requirements 

Nothing  to  report. 

Appendices 

Accepted  abstract  from  science  meeting  presentation:  EMBC15_mTBI_paper_final_rla.pdf 
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Brain  Activation  Profiles  in  mTBI:  Evidence  from  Combined 
Resting-State  EEG  and  MEG  Activity 

Lianyang  Li,  Mattia  F.  Pagnotta,  Xianghong  Arakaki,  Thao  Tran,  David  Strickland, 
Michael  Harrington,  and  George  Zouridakis,  Senior  Member,  IEEE 


Abstract — In  this  study,  we  compared  the  brain  activation 
profiles  obtained  from  resting  state  Electroencephalographic 
(EEG)  and  Magnetoencephalographic  (MEG)  activity  in  six 
mild  traumatic  brain  injury  (mTBI)  patients  and  five 
orthopedic  controls,  using  power  spectral  density  (PSD) 
analysis.  We  first  estimated  intracranial  dipolar  EEG/MEG 
sources  on  a  dense  grid  on  the  cortical  surface  and  then 
projected  these  sources  on  a  standardized  atlas  with  68  regions 
of  interest  (ROIs).  Averaging  the  PSD  values  of  all  sources  in 
each  ROI  across  all  control  subjects  resulted  in  a  normative 
database  that  was  used  to  convert  the  PSD  values  of  mTBI 
patients  into  z-scores  in  eight  distinct  frequency  bands.  We 
found  that  mTBI  patients  exhibited  statistically  significant 
overactivation  in  the  delta,  theta,  and  low  alpha  bands. 
Additionally,  the  MEG  modality  seemed  to  better  characterize 
the  group  of  individual  subjects.  These  findings  suggest  that 
resting-state  EEG/MEG  activation  maps  may  be  used  as 
specific  biomarkers  that  can  help  with  the  diagnosis  of  and 
assess  the  efficacy  of  intervention  in  mTBI  patients. 

I.  Introduction 

Mild  traumatic  brain  injury  (mTBI)  is  defined  as  a 
transient  change  in  brain  function  affecting  the  mental  state 
and  possibly  consciousness  of  patients  due  to  an  external 
force  (Menon  et  ah,  2010).  mTBI  is  difficult  to  diagnose 
because  patients  typically  lack  apparent  external  injuries  and 
clear  pathological  findings  in  conventional  computed 
tomography  and  magnetic  resonance  imaging  (MRI)  scans 
(Tarapore  et  ah,  2013),  although  evidence  of  microscopic 
MRI-based  morphological  changes  has  been  recently 
reported  (Pasternak  et  ah,  2014;  Sasaki  et  ah,  2014). 
Symptoms,  such  as  headaches,  fatigue,  and  dizziness 
(Cassidy  et  ah,  2004)  usually  emerge  on  the  day  of  injury 
and  persist  for  a  few  days  following  injury  (Boccaletti  et  ah, 
2006),  but  in  most  patients,  symptoms  resolve  and  cognition 
recovers  within  three  months.  However,  up  to  25%  of 
patients  (Sigurdardottir  et  ah,  2009)  suffer  residual 
symptoms,  long-term  impairment,  and  sometimes  disability 
(Levin  2009).  Traumatic  brain  injury  is  a  major  cause  of 
sustained  morbidity  and  disability  both  in  the  military  and 

M.  Pagnotta  and  h.  hi  are  with  the  Biomedical  Imaging  Lab, 
University  of  Houston,  Houston  TX,  USA. 

X.  Arakaki,  T.  Tran,  D.  Strickland,  and  M.  Harrington 
(mghworks@hmri.org)  are  with  Huntington  Medical  Research 
Institutes,  Pasadena,  CA. 

G.  Zouridakis  (zouridakis@uh.edu;  phone:  +1-713-743-8656; 
FAX:  +1-713-743-0172)  is  with  the  Departments  of  Engineering 
Technology,  Computer  Science,  and  Electrical  &  Computer 
Engineering,  University  of  Houston,  Houston,  TX  77204,  USA. 


civilian  populations  (Tarapore  et  ah,  2013),  as  repeated 
mTBI  can  cause  a  wide  range  of  neurological  and  cognitive 
deficits  affecting  memory,  reasoning,  language,  and 
emotions  (NINDS,  2002). 

In  the  last  two  decades,  a  plethora  of  studies  has 
attempted  to  characterize  the  structural  and  functional  effects 
of  mTBI  (Eierud  et  ah,  2014).  Of  particular  interest  are 
studies  of  resting  state  neurophysiological  recordings, 
obtained  using  Electroencephalography  (EEG)  and 
Magnetoencephalography  (MEG),  because  they  require  no 
training  or  experience  with  cognitive  tasks,  and  they  impose 
minimal  demands  on  a  patient,  which  is  especially  important 
after  brain  injury.  One  of  the  very  first  studies  using  resting- 
state  MEG  as  a  possible  biomarker  for  mTBI  showed  that 
functional  connectivity  could  be  a  valuable  tool  for  early 
detection  of  mTBI  (Zouridakis  et  ah,  2012).  Other  studies 
showed  abnormal  slowing  in  brain  areas  affected  by  TBI 
(Huang  et  ah,  2014)  and  reduced  overall  functional 
connectivity  in  TBI  patients  compared  to  controls  (Tarapore 

et  ah,  2013).  In  particular,  resting-state  MEG  source  imaging-^ - 

(Huang  et  ah,  2012)  was  able  to  detect  abnormalities  in  mild 
and  moderate  TBI  with  87%  and  100%  accuracy, 
respectively.  Furthermore,  combining  MRI  with  MEG 
(Lewine  et  ah,  1999)  could  discriminate  between  healthy 
adults  and  individuals  with  resolved  mTBI.  Compared  to 
healthy  controls,  mTBI  subjects  showed  reduced  complexity 
in  multiple  brain  areas  (Luo  et  ah,  2013).  Decreased 
connectivity  in  resting-state  MEG  may  persist  for  years  after 
mTBI  (Castellanos  et  ah,  2011),  but  the  abnormally  reduced 
connectivity  might  improve  over  time  (Tarapore  et  ah, 

2013). 

Continuing  our  earlier  attempts  to  understand  how  mTBI 
affects  communication  networks  in  the  human  brain 
(Zouridakis  et  ah,  2012,  Pollonini  et  ah,  2010,  Dimitriadis  et 
ah,  2015),  in  this  study,  we  employ  recordings  of  resting- 
state  EEG  and  MEG  and  power  spectrum  analysis  at  the 
source  level  to  investigate  abnormalities  in  brain  activation 
profiles  of  mTBI  patients. 

II.  Materials  and  Methods 

A.  Subjects 

Six  mTBI  subjects  (4  male,  2  female,  average  age  28.3±7.3) 
and  five  orthopedic  controls  (3  male,  2  female,  average  age 
29.4+7.4),  i.e.,  subjects  with  minor  orthopedic/extremity 
injuries  who  did  not  sustain  head  injury  were  recruited  for 
this  study.  All  data  were  obtained  at  the  Huntington  Medical 
Research  Institutes,  Pasadena,  CA,  USA.  Exclusion  criteria 
included  a  personal  history  of  neurological  or  psychiatric 
illness,  neurological  disorders,  serious  medical  condition. 
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and  drug  or  alcohol  addiction.  The  protocol  was  approved  by 
the  appropriate  institutional  review  board,  and  written 
informed  consent  was  obtained  from  all  participants  in  the 
study. 


B.  EEG  and  MEG  Recordings  and  Signal  Preprocessing 

Subjects  were  asked  to  remain  as  still  as  possible  during 
the  recording  procedure  and  keep  their  eyes  open. 

Approximately  5  minutes  of  continuous  EEG  activity 
was  acquired  from  each  subject  using  a  dry  electrode  EEG 
system  (Wearable  Sensing,  San  Diego,  CA).  The  system 
includes  21  channels  for  EEG  and  two  additional  channels  to 
record  EOG  and  EKG  activity  to  monitor  eye  and  cardiac 
artifacts.  Data  were  collected  at  a  sampling  rate  of  300  Hz. 

Furthermore,  5  minutes  of  continuous  MEG  activity  was 
acquired  from  each  subject  using  a  CTF  MEG  system  (MEG 
International  Services  Ltd.,  Coquitlam,  BC,  Canada).  This 
whole-head  system  includes  66  axial  gradiometer  sensors 
and  31  additional  channels  that  can  be  used  for  noise 
reduction.  Data  were  collected  at  a  sampling  rate  of  625  Hz 
and  bandpass  filtered  between  0.1-200  Hz  using  hardware. 

All  data  analyses  were  done  in  MATLAB  (The 
MathWorks,  Inc.,  Natick,  Massachusetts,  United  States) 
using  Brainstorm  (Tadel  et  ah,  2011)  and  in-house  software 
developed  in  the  Biomedical  Imaging  Lab  at  the  University 
of  Houston.  Both  the  EEG  and  MEG  sensors  were  co¬ 
registered  with  an  anatomical  MRI  template  (Colin27)  using 
six  fiducial  references,  the  nasion,  left  pre-auricular  point, 
right  pre-auricular  point,  anterior  commissure,  posterior 
commissure,  and  an  interhemisphcric  point.  All  results  were 
visually  inspected  for  accuracy. 

A  detrending  procedure  was  applied  to  all  EEG  and  MEG 
recordings  to  remove  linear  trends.  Additionally,  data  were 
bandpass  filtered  between  0.1-80  Hz,  whereas  line  noise  was 
removed  using  a  50  Hz  notch  filter.  Eye  blink  artifacts  were 
removed  using  the  additional  eye  channel  signals  and  an 
automatic  eye  blink  detection  procedure  based  on  signal- 
space  projection.  The  identified  eye  activity  topographies 
were  visually  inspected  for  accuracy.  At  least  2  minutes  of 
artifact-free  EEG/MEG  activity  was  necessary  for  a  subject 
to  be  included  in  final  analysis. 


(Hamalainen  et  ah,  1994)  and  related  software1  (MNE 
manual,  section  6).  The  dSPM  value  at  each  location  is  equal 
to  the  wMNE  value  divided  by  the  projection  of  the 
estimated  noise  covariance  matrix  onto  each  source  point. 
After  whitening,  the  operational  noise  covariance  matrix  is 
by  definition  the  identity  matrix,  and  hence  the  projection  of 
the  noise  is  equal  to  the  L2  norm  of  the  row  vector  of  the 
wMNE  inverse  operator  (in  the  case  of  fixed  dipole 
orientations). 

More  specifically,  given  a  set  of  EEG/MEG  surface 
recordings  x(t),  the  relationship  between  the  intracranial 
dipole  sources  s(t)  and  the  EEG/MEG  data  x(t)  is  given  by 
the  so  called  foi~ward  solution, 


x(t)  =  A  s(t)  +  n(t), 


(1) 


where  s(t)  denotes  a  vector  of  dipole  component  strengths,  A 
denotes  the  linear  forward  matrix  operator,  and  n(t)  denotes 
additive  noise  (Dale  and  Sereno,  1993).  Assuming  that  the 
prior  information  about  dipole  strength  follows  a 
multivariate  Gaussian  distribution,  the  maximum  a  posteriori 
probability  estimate  is  given  by 


s(t)  =  Wx{t),  with  W  =  RAt(ARAt  +  C)-1. 


(2) 


W  is  the  inverse  operator,  C  =<  n(t)n(t)r  >  is  the  data 
noise  covariance  matrix,  and  R  is  the  spatial  covariance 
matrix  of  the  dipole  strength  vector.  The  variance  of  each 
dipole  strength  estimate  due  to  the  additive  noise  n(t)  is 
given  by 


rtar(Sj)  =<  (Wj7i(t))2  >=  WjCwf, 


(3) 


where  sL  denotes  the  i  element  of  the  dipole  strength  vector 
s  and  W;  is  the  ;th  row  of  the  linear  inverse  operator  W. 

In  the  case  of  fixed  dipole  orientations,  for  each  time  point 
t  and  location  t,  a  noise-normalized  activity  estimate  sn  i(t ) 
can  be  computed  by  dividing  the  total  dipole  strength 
estimate  in  location  i  by  the  predicted  standard  error  of  the 
estimate  due  to  the  additive  noise,  using  the  formula 


C.  Intracranial  Source  Power  Analysis 

A  BEM  (boundary  element  method)  head  model  was 
used  for  the  EEG/MEG  forward  calculation.  Three  layers  of 
tissues  representing  the  scalp,  outer  skull,  and  inner  skull 
were  extracted  from  a  reference  MRI  (Colin27),  using  1922 
vertices  per  layer  and  4  mm  of  skull  thickness.  The  grid  of 
dipolar  sources  was  defined  on  the  reference  MRI  cortex 
surface  and  consisted  of  approximately  15,000  vertices.  The 
lead  field  matrix  was  computed  using  the  overlapping 
spheres  method  (Huang  et  ah,  1999). 

Estimation  of  intracranial  sources  for  each  subject  was 
performed  using  the  dynamical  Statistical  Parametric 
Mapping  (dSPM)  procedure  (Dale  et  ah,  2000),  which  is 
based  on  whitened  minimum  norm  estimation  (wMNE),  a 
depth-weighted  linear  L2  minimum  norm  estimation 
algorithm  inspired  from  the  original  MNE  method 


ViCO  = 


Siit)  wrx(t) 


y/wiCwJ1  y/wiCwJ1 


(4) 


For  source  reconstruction,  we  used  the  recordings 
obtained  from  all  sensors.  A  pre-whitening  transformation 
was  applied  to  the  data  to  pre-scale  the  channels  using  the 
noise  covariance  matrix.  Furthermore,  we  selected 
constrained  source  orientation,  which  considers  that  at  each 
vertex  of  the  cortex  surface  there  is  only  one  dipole  whose 
orientation  is  normal  to  the  surface  at  this  point. 

Power  spectral  density  (PSD)  analysis  based  on  Welch’s 
method  was  performed  on  all  source  estimates,  using  a 
window  length  of  2  sec  with  50%  time  overlap.  The 
EEG/MEG  frequency  bands  of  interest  were  the  following: 
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Fig.  3  Statistically  significant  differences  in  activation 
profiles  (z-scores)  between  an  mTBl  patient  and  the 
normative  database  in  the  lower  frequency  bands,  for  EEG 
(top)  and  MEG  (bottom). 


Fig.  1  Statistically  significant  differences  in  activation  profiles 
(z-scores)  between  mTBI  and  control  subjects  in  the  lower 
frequency  bands,  for  EEG  (top)  and  MEG  (bottom). 

0.1-4  Hz  (delta),  4-8  Hz  (theta),  8-10.5  Hz  (alphal),  and 
10.5-13  Hz  (alpha2).  For  higher  frequencies,  13-20  Hz 
(betal),  20-30  Hz  (beta2),  30-40  Hz  (gammal),  and  40-80 
Hz  (gamma2)  were  selected. 

The  Desikan-Killiany  atlas  (Desikan  et  ah,  2006)  defined 
in  the  FreeSurfer  software  was  used  for  common  co¬ 
registration  of  sources.  This  atlas  consists  of  34  brain 
regions  of  interest  (ROIs)  for  each  hemisphere  and  it  is 
available  as  a  free  download  online2.  After  averaging  the 
power  estimates  across  all  vertices  belonging  to  the  same 
ROI,  a  map  of  68  ROIs  by  8  frequency  bands  was  obtained 
for  each  subject. 

D.  Normative  Database  and  z-score  Maps 

We  further  defined  a  normative  database  using  all  five 
datasets  from  the  control  group,  separately  for  the  EEG  and 
MEG  recordings.  Specifically,  averaging  the  power  maps 
across  the  five  control  subjects  yielded  two  matrices  (68  x  8, 
ROIs-by-frequency  bands),  one  with  the  mean  values  and  a 
second  one  with  the  standard  deviations,  for  each  frequency 
band.  The  same  estimation  procedure  was  followed  for  each 
mTBI  subject. 

Each  patient  was  compared  to  the  normative  values,  and 
for  each  ROI  assessed,  a  z  score  was  computed  using  the 
following  expression: 

Pa-Mean ?• 

Zij=  \nc  with  i=  1,2...  68,  j  =  1,  2,  ...  8,  (5) 

ij 

where  PLj  is  the  ij  element  of  the  power  map  and  Meanfj 
and  SDfj  are  the  mean  and  standard  deviation  values, 
respectively,  from  the  two  normative  database  matrices 
corresponding  to  the  i-th  ROI  and  y'-th  frequency  band 
(Huang  et  al.,  2012). 

E.  Statistical  Thresholding 

To  identify  statistically  significant  z-scores  we  used 
alpha  =  0.05  and  false  discovery  rate  correction  (Benjamini  & 

2  http://surfer.nmr.mgh.harvard.edu/ 


Fig.  2  Statistically  significant  differences  in  activation 
profiles  (z-scores)  between  a  control  subject  and  the  normative 
database  in  the  lower  frequency  bands,  for  EEG  (top)  and 
MEG  (bottom). 


Hochberg,  1995)  for  the  544  (68  x  8)  multiple  comparisons, 
resulting  in  a  threshold  value  of  2.2421.  After  thresholding, 
the  statistically  significant  z-scores  were  projected  onto  the 
Desikan-Killiany  atlas  (Desikan  et  al.,  2006)  using  different 
colors  for  each  z-score  value  to  visualize  the  brain  areas  that 
differed  significantly  between  the  two  groups. 

III.  Results 

When  comparing  the  mTBI  and  control  groups, 
statistically  significant  differences  in  the  form  of 
overactivation  were  seen  primarily  in  the  theta  (4-8  Hz)  and 
low  alpha  (8-10.5  Hz)  bands  for  the  EEG  and  the  delta  (0.1- 
4  Hz),  theta  (4-8  Hz),  and  low  alpha  (8-10.5  Hz)  bands  for 
MEG  data,  as  shown  in  Fig.  1. 

To  test  the  power  of  the  methodology  to  correctly 
classify  subjects  on  a  single  subject  basis,  we  used  the  first 
four  controls  to  construct  the  normative  database,  while  the 
fifth  one  was  used  to  construct  the  activation  maps  (z  scores) 
shown  in  Fig.  2.  For  comparison,  Fig.  3  shows  the  activation 
maps  of  one  mTBI  subject.  These  figures  show  that  the 
control  subject  does  not  differ  significantly  from  the 
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normative  database,  while  the  mTBI  subject  shows 
significant  differences  mostly  in  the  lower  frequency  bands. 
Additionally,  the  MEG  modality  seems  to  more  accurately 
characterize  the  individual  subject  group. 

Overall,  even  though  the  number  of  participants  is  small 
at  this  early  stage  of  the  study,  the  results  obtained  suggest 
that  analysis  of  resting-state  EEG  and  MEG  activation  maps 
is  a  powerful  tool  that  can  help  in  the  diagnosis  of  and  assess 
the  efficacy  of  intervention  in  mTBI. 

Acknowledgment 

This  work  was  supported  in  part  by  a  grant  from  the 
Department  of  Defense  to  HMRI,  U.S.  Army  W81XWH-13- 
1-0005. 

References 

[1]  Benjamini,  Y.,  &  Hochberg,  Y.  (1995).  Controlling  the  false 
discovery  rate:  a  practical  and  powerful  approach  to  multiple 
testing.  Journal  of  the  Royal  Statistical  Society.  Series  B 
(Methodological),  289-300. 

[2]  Boccaletti,  S.,  Latora,  V.,  Moreno,  Y.,  Chavez,  M.,  &  Hwang, 
D.  U.  (2006).  Complex  networks:  Structure  and  dynamics. 
Physics  reports,  424(4),  175-308. 

[3]  Cassidy,  J.  D.,  Carroll,  L.,  Peloso,  P.,  Borg,  J.,  Von  Holst,  H., 
Holm,  L.,  ...  &  Coronado,  V.  (2004).  Incidence,  risk  factors 
and  prevention  of  mild  traumatic  brain  injury:  results  of  the 
WHO  Collaborating  Centre  Task  Force  on  Mild  Traumatic! 
Brain  Injury.  Journal  of  rehabilitation  medicine,  36(0),  28-60. 

[4]  Castellanos,  N.P.  (2011).  Principles  of  recovery  from 
traumatic  brain  injury:  Reorganization  of  functional  networks. 
Neuroimage,  55,  1189-1199. 

[5]  Dale,  A.  M.,  Liu,  A.  K„  Fischl,  B.  R.,  Buckner,  R.  L., 
Belliveau,  J.  W.,  Lewine,  J.  D.,  &  Halgren,  E.  (2000). 
Dynamic  statistical  parametric  mapping:  combining  fMRI  and 
MEG  for  high-resolution  imaging  of  cortical 
activity.  Neuron,  26(1),  55-67. 

[6]  Dale,  A.  M.,  Sereno,  M.  I.  (1993).  Improved  localization  of 
cortical  activity  by  combining  EEG  and  MEG  with  MRI 
cortical  surface  reconstruction:  a  linear  approach.  Journal  of 
Cognitive  Neuroscience,  5,  162-176. 

[7]  Desikan,  R.  S.,  Segonne,  F.,  Fischl,  B.,  Quinn,  B.  T., 
Dickerson,  B.  C.,  D.  Blacker,  Buckner,  R.  L.,  Dale,  A.  M., 
Maguire,  R.  P.,  Hyman,  B.  T.,  Albert,  M.  S.,  Killiany,  R.  J. 
(2006).  An  automated  labeling  system  for  subdividing  the 
human  cerebral  cortex  on  MRI  scans  into  gyral  based  regions 
of  interest.  Neuroimage,  31,  968-980. 

[8]  Dimitriadis,  S.,  Zouridakis,  G.,  Rezaie,  R.,  Babajani-Feremid, 
A.,  and  Papanicolaou,  A.  C.  (2015).  Functional  Connectivity 
Changes  Detected  with  Magnetoencephalography  after  Mild 
Traumatic  Brain  Injury.  Neuroimage:  Clinical  (under  review). 

[9]  Eierud,  C.,  Craddock,  R.C.,  Fletcher,  S.,  Aulakh,  M.,  King- 
Casas,  B.,  Kuehl,  D.,  LaConte,  S.M.  (2014).  Neuroimaging 
after  mild  traumatic  brain  injury:  Review  and  meta-analysis. 
Neuroimage:  Clinical,  4,  283-294. 

[10]  Engemann,  D.  A.,  Gramfort,  A.  (2015).  Automated  model 
selection  in  covariance  estimation  and  spatial  whitening  of 
MEG  and  EEG  signals.  Neuroimage,  108,  328-342. 

[11]  Hamalainen,  M.  S.,  &  Ilmoniemi,  R.  J.  (1994).  Interpreting 
magnetic  fields  of  the  brain:  minimum  norm 
estimates.  Medical  &  biological  engineering  &  computing, 
32(1),  35-42. 

[12]  Huang,  M.  X.,  Mosher,  J.  C.,  Leahy,  R.  M.  (1999).  A  sensor- 
weighted  overlapping-sphere  head  model  and  exhaustive  head 


model  comparison  for  MEG.  Physics  in  Medicine  and 
Biology,  44(2),  423-440. 

[13]  Huang,  M.  X.,  Nichols,  S.,  Baker,  D.  G.,  Robb,  A.,  Angeles, 
A.,  Yurgil,  K.  A.,  ...  &  Lee,  R.  R.  (2014).  Single-subject-based 
whole-brain  MEG  slow-wave  imaging  approach  for  detecting 
abnormality  in  patients  with  mild  traumatic  brain  injury. 
Neuroimage:  Clinical,  5,  109-119. 

[14]  Huang,  M.  X.,  Nichols,  S.,  Robb,  A.,  Angeles,  A.,  Drake,  A., 
Holland,  M.,  Asmussen,  S.,  ...  &  Lee,  R.  R.  (2012).  An 
automatic  MEG  low-frequency  source  imaging  approach  for 
detecting  injuries  in  mild  and  moderate  TBI  patients  with  blast 
and  non-blast  causes.  Neuroimage,  61(4),  1067-1082. 

[15]  Levin,  H.S.  (2009).  Mission  Connect  Mild  TBI  Translational 
Research  Consortium.  Baylor  College  of  Medicine  Houston 
TX. 

[16]  Lewine,  J.D.,  Davis,  J.T.,  Sloan,  J.H.,  Kodituwakku,  P.W., 
Orrison,  W.W.  (1999).  Neuromagnetic  assessment  of 
pathophysiologic  brain  activity  induced  by  minor  head  trauma. 
American  Journal  of  Neuroradiology,  20(5),  857-866. 

[17]  Luo,  Q.,  Xu,  D.,  Roskos,  T.,  Stout,  J.,  Kull,  L.,  Cheng,  X.,  ... 

&  Bucholz,  R.  D.  (2013).  Complexity  analysis  of  resting  state 
magnetoencephalography  activity  in  traumatic  brain  injury 
patients.  Journal  of  neurotrauma,  30(20),  1702-1709. 

[18]  Menon,  D.  K.,  Schwab,  K„  Wright,  D.  W„  &  Maas,  A.  I. 
(2010).  Position  statement:  definition  of  traumatic  brain 
injury.  Archives  of  physical  medicine  and  rehabilitation, 
91(11),  1637-1640. 

[19]  N1NDS  National  Institute  of  Neurological  Disorders  and 
Stroke.  Traumatic  brain  injury:  hope  through  research. 
Bethesda  (MD):  National  Institutes  of  Health,  NIH  Publication 
No.  02-2478,  2002. 

[20]  Pasternak,  O.,  Koerte,  1.  K.,  Bouix,  S.,  Fredman,  E.,  Sasaki, 

T.,  Mayinger,  M.,  ...  &  Echlin,  P.  S.  (2014).  Hockey  / 
Concussion  Education  Project,  Part  2.  Microstructural  white  ^ 
matter  alterations  in  acutely  concussed  ice  hockey  players:  a 
longitudinal  free-water  MRI  study:  Clinical  article.  Journal  of 
neurosurgery,  120(4),  873-881. 

[21]  Pollonini,  L„  Pophale,  S.,  Situ,  N„  Wu,  M.  H.,  Frye,  R.  E, 
Leon-Carrion,  J.,  &  Zouridakis,  G.  (2010).  Information 
communication  networks  in  severe  traumatic  brain 
injury.  Brain  topography,  23(2),  221-226. 

[22]  Sasaki,  T.,  Pasternak,  O.,  Mayinger,  M.,  Muehlmann,  M., 
Savadjiev,  P.,  Bouix,  S.,  ...  &  Koerte,  I.  K.  (2014).  Hockey 
Concussion  Education  Project,  Part  3.  White  matter 
microstructure  in  ice  hockey  players  with  a  history  of 
concussion:  a  diffusion  tensor  imaging  study:  Clinical 
article.  Journal  of  neurosurgery,  120(4),  882-890. 

[23]  Sigurdardottir,  S.,  Andelic,  N.,  Roe,  C.,  Jerstad,  T.,  & 
Schanke,  A.  K.  (2009).  Post-concussion  symptoms  after 
traumatic  brain  injury  at  3  and  12  months  post-injury:  a 
prospective  study.  Brain  Injury,  23(6),  489-497. 

[24]  Tadel,  F.,  Baillet,  S.,  Mosher,  J.  C.,  Pantazis,  D.,  Leahy,  R.  M. 
(2011).  Brainstorm:  A  User-Friendly  Application  for 
MEG/EEG  Analysis.  Computational  Intelligence  and 
Neuroscience,  2011:  879716. 

[25]  Tarapore,  P.  E.,  Findlay,  A.  M.,  LaHue,  S.  C.,  Lee,  H., 
Honma,  S.  M.,  Mizuiri,  D.,  ...  &  Mukherjee,  P.  (2013). 
Resting  state  magnetoencephalography  functional  connectivity 
in  traumatic  brain  injury.  Journal  of  neurosurgery,  1 18(6), 
1306. 

[26]  Zouridakis,  G.,  Patidar,  U.,  Situ,  N.,  Rezaie,  R.,  Castillo,  E. 
M„  Levin,  H.  S.,  &  Papanicolaou,  A.  C.  (2012).  Functional 
connectivity  changes  in  mild  traumatic  brain  injury  assessed 
using  magnetoencephalography.  Journal  of  Mechanics  in 
Medicine  and  Biology,  12(02). 


